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The  Angular  Impact  Distibution  of 
Charged  Particles  Attracted  to  a 
Charged  Cylindrical  Spacecraft 


1.  INTRODUCTION 

Spacecraft  in  geosynchronous  orbit  exhibit  charging  when  immersed  in  hot 
substorm  plasmas.  The  potential  that  a  surf-ice  attains  depends  both  on  the 
material's  properties  and  on  the  plasma's  characteristics.  Secondary  emission 
of  electrons  by  electrons  and  ions,  as  well  as  backsrattering,  determine  the 
floating  potential.  These  are  both  functions  of  the  angle  of  incidence  of  particles 
on  the  surface. 

In  this  work,  the  angular  distribution  of  incident  particles  is  obtained  for 
attracted  and  repelled  particles  for  cylindrical  spacecraft.  Results  are  given 
for  the  plane  and  spherical  cases  as  well. 

For  isotropic  spacecraft  materials,  the  backscattering  coefficients,  the 
secondary  emission  coefficients,  and  the  shielding  efficiency  depend  on  impact 
energy  and  impact  zenith  angle,  but  are  independent  of  azimuth  angle.  According¬ 
ly,  only  the  distribution  of  impact  zenith  angles  is  important.  For  convenience, 
only  protons  will  be  considered,  the  extension  to  other  charged  particles  being 
obvious.  First  the  problem  will  be  precisely  defined,  then  a  simplified  model  of 
physical  reality  set  up,  and,  finally,  the  distribution  of  impact  angles  derived. 
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2.  DF.FIMTION  OF  PKOlil.FM 


Kind  the  flux,  j  of  protons  striking  an  element  of  ami  £A  with  impart 
kinrtir  energies  within  tho  band  K  t  J1K2  and  with  impart  zrnith  angles  witlitn 
tho  band  0  t  £0/2. 


X  MODF.I. 


The  model  is  that  of  an  infinitivelv  long  rvlindor  of  radius  H  .  Beyond  a 

l  ) 

"sheath"  of  radius  11^,  tho  plasma  is  isotropio  and  is  at  zero  potential.  1'ho 
potential  within  tho  sheath  is  modeled  bv: 


ln(K/R  .) 

o  ln(R  7R  ) 

o  s 


(1) 


The  subscripts  "o"  and  "s"  will  bo  used  to  indicate  values  of  variables  at  H  and 

R  No  subscript  will  be  used  for  general  values  of  R. 

&  -a  -1-1-1 

In  tho  model,  the  differential  flux  f(K,  protons  cm  “  see  ster  eV  is 

a  function  of  kinetic  and  potential  energies  only  —  except  in  the  "shadow"  of  the 

Satellite  where  it  is  zero.  Both  potential  ($)  and  kinetic  (K)  energy  are  expressed 

in  electron  volts.  Any  impact  angle  and  energy  corresponding  to  a  trajectory 

originating  elsewhere  on  the  satellite  is  considered  to  bo  "shadowed.” 


4.  DERIVATION  OF  IMPACT  /.F.NITII  ANGLE  DISTRIIUTION 

The  impact  proton  flux  is  the  product  of  four  factors,  namely:  energy  band 
width,  unshadowed  solid  angle,  projected  area,  and  differential  flux.  This  may 
be  w  ritten  as 

j  -  {£i'.)(£0  4  sin  0  )(£A  cos  0  )[f( K  ,  $  ))  (2) 

•  l>  7.  O  O  O'  o 

where  4^  =  total  unshadowed  azimuth  angle.  The  solid  angle  is  that  which  would 

be  subtended  at  the  earth's  center  by  the  area  enclosed  by  latitudes  0  ±  £0:2  and 

'  "  o 

longitudes  4  i  4z/2.  The  azimuth  angle  will  be  taken  as  zero  perpendicular  to  the 
cylinder's  axis.  If  shadowing  occurs,  it  will  be  at  the  larger  angles. 

Equation  (2)  may  be  rewritten  as 

j  fAE  £0  £A  f(E  ,  *  )l[sin  0  cos  9)144  I  (3) 

Jo  o  o  o  o  ' 
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where  -  maximum  unshadowed  azimuth  angle,  <t>  being  measured  So  as  always 
to  lie  in  the  first  quadrant.  The  factor  4  accounts  for  the  four  quadrants.  A 
further  rewriting  yields 

i  =  K  <f>  sin  2  0 
■o  tn  o 

K  =  ^  AO  AA  f(E  ,  $  )  (4) 

2  <>  o 

0  <  0  <  n  f  2  . 

""  m 

Note:  sin  2  0  =2  sin  0  c os  0 
o  o  o 

Equation  (4)  is  very  general.  The  geometry  of  the  spacecraft's  surface  and  field 
enter  into  the  equation  only  through  <J|  ,  the  maximum  unshadowed  azimuth  angle. 

5.  SPECIAL  CASES 

For  positive  spacecraft  potentials,  the  reverse  trajectories  curve  away  from 
the  spacecraft;  therefore,  there  is  no  shadowing  and 

6  ---  *12  if  $  ?-  0  .  (a) 

m  ' 

From  Eq.  (4),  the  angular  distribution  is: 

j  =  (tr/2)  K  sin  2  0  . 

■  o  ' 

The  other  special  case  is  that  of  looking  for  protons  of  impossibly  low  energy. 
Here  we  may  write 

<j>  -  0  if  9  <  -E  ,  ( f>) 

m  o  o 

which  has  meaning  only  for  negatively  charged  spacecraft. 


6.  MAXIMUM  UNSHADOWED  AZIMUTH  ANGLE 

The  simplest  cases  have  been  dealt  with  in  the  preceding  section.  This 
section  will  consider  cases  defined  by: 


(7) 


-E  <  $>  <0 

o  o 


The  model  has  zero  axial  field.  Therefore,  the  kinetic  energy  associated  with 
the  axial  velocity  component  is  independent  of  radial  position  and  may  be  written: 

E  -  E  sin2  <j>  sin2  6  ,  (8) 

a  o  o  o 

where  E  is  the  impact  energy  and  <4o  and  8^  are  the  impact  angles.  This  rela¬ 
tionship  follows  from  the  breakdown  of  the  impact  velocity  into  four  components, 
as  shown  in  Figure  1.  The  model  has  a  radial  field  about  the  cylinder's  axis. 
Thus  the  tangential  (l  to  radial  and  axial  directions)  kinetic  energy  is  subject  to 
conservation  of  angular  momentum  and  is  given  by: 

E  (R)  =  E  (R  /R)2  cos2  <4  sin2  6  .  (9) 

t  o  o  o  o 

The  radial  kinetic  energy  is  simply: 


Ef  =  E  cos"  0  .  (10) 

The  total  proton  energy  is  independent  of  the  radius,  and  beyond  the  sheath  it 
equals  kinetic  energy.  The  total  energy  is  thus: 

E  =  E+  E+Ef  +  $  ,  <  0  .  (11) 

rat 

At  the  sheath's  edge  Eq.  (11)  becomes: 

E  =  E  cos2  0  +  E  sin2  <4  sin2  0  +  (R  /R  )2  E  cos2  <4  sin2  0  (12) 

s  so  o  o  o'  S  o  o  o 

and  at  the  satellite's  surface  it  becomes 

E  =  E  cos2  0  +  E  sin2  $  sin2  8  +  E  cos2  <4  sin2  9  +  $  (13a! 

o  oo  o  oo  o  oo 


E  =  E  +  $ 
s  o 


$  <  0  . 


The  limiting  case  corresponds  to  a  vanishing  radial  energy  at  the  sheath's  edge  or 


cos  0=0 
s 


I 
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Equations  (4)  and  (li>)  together  give  the  zenith  angle  (0)  distribution  of  protons 
striking  the  surface  of  a  cylindrical  satellite. 

The  angular  distribution  derivation  given  has  considered  only  the  potential  at 
the  inner  (R  -  R  )  and  outer  (R  =  Rg)  edges,  not  the  potential  within  the  sheath 
(R  <  R  ,).  It  may  be  shown  that  the  potential  within  the  sheath  is  noi  controlling, 

s  '  i 

provided  that  at  all  points  within  the  sheath  it  is  equal  to  or  more  negative  than 
that  of  the  inverse  square  potential  function: 


<t> 


(R  /R)2- 1 
s _ 

( r  / r  r- i 


<  0 

o 


(17) 


This  function  is  adjusted  to  match  the  model  at  the  sheath's  edge.  The  logarith¬ 
mic  function  of  the  model  satisfies  this  condition.  The  proof  is  similar  to  that 
given  for  the  spherical  case  for  objects  in  a  gravitational  well.  The  proof  is  only 
valid  if  the  potential  is  a  function  of  radial  position  alone. 

The  distribution  of  impact  zenith  angles  according  to  Eqs.  (4)  and  (lii)  is 
shown  in  Figure  2  for  several  cases.  The  deviations  from  the  symmetrical 
sin  (2  0)  distribution  are  due  entirely  to  shadowing.  The  shadowing  is  equivalent 
to  a  "horizon"  at  less  than  90°  below  the  zenith  due  to  refraction  of  the  protons. 


0°  10°  20°  30°  40®  50®  60°  70°  80®  90° 


ZENITH  ANGLE 

Figure  2.  Distribution  of  Zenith  Impact  Angles  for  Protons  on  a  Cylindrical 
Satellite:  (a)  Thin-Sheath  C  ase 
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T.  DISCI  SS|()N 


It:,-  •hf.iy  given  here  is  cx:iil  in  !  1 1  n  t  n  is  completely  n  .nsis'  •  iili  the 

model.  The  mi >del,  hmvrvci',  i epresonts  a  physical  impossibility  in  that  tin  * -1  •*»  - 
*  r  i  e  field  is  discont  i minus  ;it  fin’  surface  of  the  sin-nth  and  an  approximation  in 
that  inici  effects  are  taken  to  be  vanishingly  small.  Only  if  the  nearest  end  is  at  a 
distance  very  large  compared  to  the  radius  will  its  effect  heroine  vanishingly 
small.  file  theory  does,  however,  represent  a  limiting  ease  and,  as  such,  is 
useful. 


8.  SIMMVIO 


The  angular  distribution  of  impacting  repelled  particles  is  sin  (2  »),  as  given 
by  Kq.  (4),  since  there  is  no  shadowing  for  repelled  particles,  electrons  inci¬ 
dent  .-n  a  negatively  charged  surface  and  ions  on  a  positively  charged  surface  w  ;l) 
impact  in  this  way.  In  order  to  calculate  surface  potentials,  taking  account  -f 
secondary  electron  emission,  one  must  know  the  distribution  of  impact  angles  aB 
well  : : s  the  secondary  electron  yield  as  a  function  of  angle. 

The  angular  distribution  of  impacting  attracted  particles  (electrons  on  a 
positive  surface,  and  positive  ions  on  a  negative  surface)  is  more  complicated. 

The  m parting  angular  distribution  for-  a  cylindrical  surface  for  attracted  parte  1,  -■ 
has  been  derived  here,  as  shown  in  Figures  2  (a),  2  (b),  and  2  (e),  for  a  range-  of 
energies  and  sheath  thickness. 

The  distribution  of  impact  angles  (attracted  particles)  for  a  plane  surface  is 

'll*  same  as  'he  case  of  zero  sheath  thickness.  For  this  rase  the  distribution  is 

proportional  to  sin  (2  0)  for  0  less  than  0  ,  where 

m 

-  cos'1  (4>/K)1/2  , 

and  where  K  is  the  impact  energy.  If  6  is  greater  than  0  the  distribution  func¬ 
tion  is  zero.  For  spherical  surfaces  there  is  no  shadowing,  so  the  distribution  :s 
again  proportional  to  sin  2  0  in  the  thick  sheath  limit. 
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